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Abstract 

The  concept  of  a  new  test  facility  to  evaluate  the  un¬ 
derwater  launch  and  flight  performance  of  full-scale  mis¬ 
siles  and  associated  launcher  systems  is  described.  Re¬ 
sults  from  preliminary  studies  on  estimates  of  the  platform 
dimensions,  hydrodynamic  characteristics,  and  platform 
response  to  a  missile  launch  impulse  are  included.  The 
VS  AERO  panel  method  code  was  used  to  generate  hy¬ 
drodynamic  force  and  moment  coefficients  for  the  condi¬ 
tion  of  steady  tow  at  5  knots.  These  coefficients  were  then 
used  to  evaluate  the  platform  trim  stability  under  steady 
tow  conditions.  Calculations  indicate  that  the  platform  is 
stable  during  all  three  phases  of  operation:  missile  in  tube, 
empty  tube,  and  tube  fully  flooded.  Preliminary  estimates 
of  the  platform  response  due  to  a  launch  impulse  show  a 
recoil  that  is  not  expected  to  significantly  alter  the  fidel¬ 
ity  of  the  in-tube  trajectory  data. 

1.0  Introduction 

Since  the  1950s,  the  U.S.  Navy  has  successfully  de¬ 
veloped  a  number  of  strategic  and  tactical  missile  weapon 
systems  that  are  launched  vertically  from  submerged,  for¬ 
ward-moving  submarines.  These  weapon  systems  cover 
a  wide  spectrum  of  designs.  These  designs  range  from 
the  Trident  II  intercontinental  ballistic  missile  to  the  Toma¬ 
hawk  Vertical  Launch  System  (VLS).  The  Trident  II  mis¬ 
sile  is  nearly  7  feet  in  diameter,  weighs  over  100,000 
pounds,  and  is  launched  from  a  submarine  moving  at  near¬ 
zero  forward  speed,  while  the  Tomahawk  VLS,  which  is 
21  inches  in  diameter  in  torpedo  tube,  weighs  less  than 
4000  pounds  and  is  capable  of  being  launched  at  moder¬ 
ate  forward- way  with  maneuvering. 
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The  successful  development  of  these  complex  weapon 
systems  is  in  no  small  measure  a  result  of  the  ability  to 
test,  evaluate,  and  understand  the  physics  of  underwater 
launch.  Historically,  the  engineering  development  of  a  new 
submarine-launched  missile  combines  both  computer 
simulation  and  subscale  or  full-scale  test  with  approxi¬ 
mation  to  the  actual  launch  conditions.  While  some  launch 
phenomena  can  be  mimicked  in  this  manner,  complex 
interactions  between  the  missile,  launcher,  and  the  en¬ 
vironment  can  be  missed.  Extremely  costly  “fixes”  and 
programmatic  delays  can  arise  if  complete  weapon  sys¬ 
tem  testing  is  delayed  until  deployed  on  the  actual  sub¬ 
marine.  The  need  for  high-fidelity  and  robust  testing  dur¬ 
ing  the  developmental  phase  of  a  missile  is  important; 
however,  this  need  is  tempered  with  the  constraint  to  con¬ 
trol  costs  associated  with  developing  and  maintaining  a 
full-scale  launch  test  facility. 

One  potential  solution  to  this  ever-present  engineer¬ 
ing  dilemma  is  a  concept  called  the  Towed  Underwater 
Launch  Platform  (TULP).  The  TULP  is  a  submergible, 
towed  platform  that  simultaneously  replicates  the  hydro¬ 
static  pressure  (depth)  and  the  submarine  deck  flow  field 
(crossflow  velocity)  at  full  scale  to  capture  the  complex 
phenomenology  of  underwater  launch.  Sufficient  flexibil¬ 
ity  is  built  into  the  design  to  accommodate  a  diverse  range 
of  missile  types  and  test  requirements. 

2.0  TULP  Concept  Description 

The  overall  guidance  in  developing  the  TULP  con¬ 
cept  was  finding  a  low-cost,  full-scale  facility  for  mis¬ 
siles  and  launchers  undergoing  engineering  development. 
References  1-3  document  the  initial  study  results  that 
defined  the  TULP  concept  and  explored  a  point  design 
from  the  perspective  of  physical  features  (weight,  mo¬ 
ments  of  inertia,  displacement,  and  reserve  buoyancy), 
hydrodynamic  characteristics  (drag,  tow  horsepower 
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requirement,  impulse  response,  etc.),  a  preliminary  depth 
control/deballasting  system,  and  an  operational  scenario. 

The  TULP  configuration  is  shown  in  Figures  1  and 
2.  The  principal  TULP  components  consist  of  a  support 
frame  with  deck  plating,  a  launch  tube  housed  within  an 
outer  streamlined  structure,  a  multiple  pontoon  buoyancy 
system,  and  a  closed-loop  hovering/deballasting  system. 
In  addition,  four  free-flooding  tubular  supports  between 
the  deck  frame  and  the  mount  tube  are  shown  to  provide 
lateral  and  torsional  mount  tube  rigidity.  The  rectangular 
form  of  the  deck  has  high  hydrodynamic  inertia  (added 
mass)  in  the  vertical  plane  to  minimize  platform  recoil  to 
the  launch  pulse.  The  structure  was  sized  to  control  plat¬ 
form  recoil  from  the  launch  pulse  with  limitation  to  allow 
the  system  to  operate  from  existing  Navy  wharf  facilities. 

TULP  buoyancy  is  provided  by  a  hard-tank  pontoon 
system  secured  to  the  platform  deck  and  by  the  faired 
compartment  surrounding  the  launch  tube.  The  two  large 
pontoons,  outboard  of  the  launch  tube,  serve  as  seawater 
blow/flood  tanks  and  as  a  compartment  to  house  valves 
and  pressure  sensors  associated  with  the  closed-loop 
depth-control/deballasting  system. 

The  free-flooding  structure  (“fairwater’O  added  to  the 
main  deck  replicates  the  flow  over  the  submarine  near 
the  launch  tube  region.  The  fairwater  deck  is  shown  asym¬ 
metric  about  the  platform  centerline.  To  help  reduce  any 
yaw  motion  that  would  be  induced  by  dynamic  pressure 
asymmetry  during  tow  operations,  a  flow-through  path  is 
shown. 

The  structure  encircling  the  launch  tube  provides  a 
watertight  compartment  to  house  the  launch  tube,  gas 
generator,  and  test  instrumentation.  The  compartment  al¬ 
lows  access  to  launch  equipment  during  surface  opera¬ 
tions.  A  free-flooding  fairing  is  added  to  the  encircling 
structure  to  minimize  drag  for  reduced  tow  power  during 
crossflow  tests.  Access  to  the  housing  is  provided  through 
a  service  hatch  located  under  the  fairwater  deck.  The  lower 
portion  of  the  faired  structure  contains  lead  ballast  to  aid 
in  platform  stability.  Forward- way  is  provided  by  means 
of  a  support  Navy  tugship  that  is  outfitted  with  launcher 
and  fire  control  subsystems.  A  tow  cable  with  strength 
member,  electrical  power,  and  communication  line  would 
provide  a  safe  standoff  distance  between  the  TULP  and 
towship.  In  keeping  with  past  safe  practice,  an  inert  launch 
test  vehicle  (LTV)  would  be  the  test  missile  employed  for 
engineering  development  of  any  new  system. 

Table  1  is  a  compilation  of  the  results  of  a  recent 
TULP  sizing  study.  Preliminary  estimates  of  the  TULP 
component  structure  were  generated  using  standard  ma¬ 
rine  design  practices  for  calm-to-moderate  sea  conditions 
(it  is  assumed  that  engineering  testing  would  not  be  per¬ 
formed  under  high  sea-state  conditions).  The  pontoon  size 
and  lead  ballast  were  selected  to  maintain  a  positive  sub¬ 
merged  metacentric  height  under  the  following  conditions: 


Table  1.  Estimated  principal  characteristics  of  Towed  Under- 


water  Launch  Platform. 

Principal  Dimensions 

Overall  length,  ft 

60.0 

Overall  beam,  ft 

40.0 

Overall  height,  ft 

47.5 

Surfaced  draft,  ft 

39.5 

Physical  Characteristics 

Weight  (excluding  test  vehicle),  lb 

360,000 

Submerged  displacement  (excluding 

free  flooding),  lb 

490  ,000 

Reserve  buoyancy,  lb 

24,000 

Roll  and  pitch  moment  of  inertia,  slug-ft^ 

3.0  X  10^ 

Center  of  gravity  (C.G.),*  ft 

-8.8 

Hydrodynamic  Characteristics 

Center  of  buoyancy  (C.B.),*  ft 

-4.2 

Submerged  metacentric  height  (C.B  -  C.G.),  ft 

4.6 

Forward-way  drag  @  5  knots,  lb 

-6,500 

Submerged  tow  @  5  knots,  hp 

-100 

Heave  added  mass  (including  free  flooding). 

slugs 

1.65  X  10^ 

Roll  and  pitch  added  moment  of  inertia, 

slug-ft^ 

6.25  X  10^ 

Roll  and  pitch  natural  period,  s 

18  to  19 

Depth  Control/Deballasting  System 

Ballast  tank  volume  (each),  ft^ 

850 

Ballast  tank  capacity  (each),  lb 

55,000 

Maximum  flow  rate  @  DP  =  50  PSID,  Ib/s 

540 

*The  vertical  positions  of  C.G.  and  C.B.  are  defined  with  refer¬ 
ence  to  the  axis  system  shown  in  Figure  1.  C.G.  is  with  test 
vehicle  in  tube. 


a  missile  in  the  launch  tube,  an  empty  launch  tube,  and  a 
launch  tube  filled  with  seawater. 

3.0  Hydrodynamic  Analysis 

3.1  VSAERO  Panel  Code 

This  paper  addresses  the  hydrodynamic  and  stability 
issues  associated  with  steady-state  tow.  In  order  to  assess 
the  stability  characteristics  of  the  TULP  in  steady  tow,  a 
hydrodynamic  analysis  was  performed  using  a  low-order 
panel  method  code,  VSAERO."^’^  The  code  calculates  the 
inviscid,  linearized  potential  flow  external  to  a  body.  An 
integral  boundary-layer  method  may  be  coupled  with 
VSAERO  to  obtain  an  estimate  of  the  viscous  effects.  For 
incompressible  and  irrotational  flow,  the  potential  flow 
will  satisfy  Laplace’s  equation.  The  potential  flow  must 
also  satisfy  the  integral  equations  obtained  from  the  appli¬ 
cation  of  Green’s  theorem.  VSAERO  solves  for  the  flow 
potential  by  solving  the  integral  equation  from  Green’s 
theorem  on  the  surface  boundary  using  a  constant  strength 
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distribution  of  sources  and  doublets  along  a  defined  set 
of  body  panels.  The  panel  method  solution  is  obtained  by 
solving  the  integral  equation  as  a  set  of  equations  and 
unknowns  equal  to  the  number  of  panels  defining  the  body 
surface.  For  the  solution  of  the  TULP  configuration,  ini¬ 
tially  a  zero  normal  velocity  boundary  condition  was  ap¬ 
plied  on  each  body  panel. 

The  vorticity  shed  by  the  body  into  the  flow  is  mod¬ 
eled  by  a  wake  off  the  trailing  edge  of  the  body.  The  wake 
is  also  modeled  by  a  set  of  panels,  by  defining  wake  grid 
planes  and  wake  lines.  The  wake  is  allowed  to  relax  or 
conform  to  a  steady-state  shape  that  satisfies  the  normal 
velocity  boundary  condition  across  the  wake  panels.  The 
integral  boundary-layer  method  can  compute  laminar  and 
turbulent  attached  flow,  boundary-layer  transition  and 
separation  points,  and  relaminarization.  Application  of  the 
integral  boundary-layer  method  adds  a  transpiration  term 
in  the  boundary  condition  to  account  for  the  boundary- 
layer  displacement  thickness.  The  data  were  postprocessed 
with  the  OMNI3D  code,^  and  wake  preprocessing  was 
performed  using  the  SPIN  code.^ 

The  analysis  was  performed  using  the  computer  fa¬ 
cilities  at  the  U.S.  Navy  Hydrodynamic/Hydroacoustic 
Technology  Center  at  the  Naval  Surface  Warfare  Center, 
Carderock  Division,  Bethesda,  Maryland.  The  calculations 
were  performed  on  a  Silicon  Graphics  Power  Challenger 
with  six  75-MHz  processors  and  2048  Mbytes  of  memory. 

The  hydrodynamic  coefficients  computed  by  the 
VS  AERO  code  were  used  to  evaluate  the  stability  of  the 
platform  in  steady  tow.  The  following  section  describes 
the  approach  used  to  develop  the  fully  submerged  plat¬ 
form  trim  stability  under  steady  tow. 

3.2  Trim  Equations  for  the  TULP  vehicle 

Reference  8  contains  the  equations  of  motion  for  a 
submarine.  For  the  TULP  trim  analysis,  the  components 
for  towcable  tension  and  towpoint  location  were  added. 
Also,  the  force  and  moment  coefficient  definitions  were 
replaced  with  the  nomenclature  consistent  with  VS  AERO 
nomenclature.  In  addition,  the  following  assumptions  were 
made:  steady  horizontal  flight,  all  coupled  hydrodynamic 
terms  set  to  zero,  all  acceleration  and  rotational  rates 
nulled,  and  zero  sideslip  and  roll  angles.  With  these  as¬ 
sumptions,  the  following  set  of  equations  result  for  the 
axial  and  normal  force  and  the  pitching  moment: 

Axial  force  equation 

|pM2/lQ+(W-B)a-7i=0  (1) 

Normal  force  equation 

^pu'^AiC^  +  C^„a) -  {W-B)-T,  =  0  (2) 


Pitching  moment  equation 

^pu^A{C^  +  C„„a)  +  (xoW-xbB) 

2  (3) 

-  (zgW  -  ZbB)ol  +  TjcZtp  -  T^xjp  =  0 

Equations  1,  2,  and  3  are  solved  simultaneously  to 
obtain  the  equation  for 

1  0 

—  pM  A(Cfnl-  Cj\ztp  +  Cp^xjp)  -  {xqW  ~  xpB)  +  xjp{W  -B) 
-^pu^AiC^l  -  Cho^Xtp)  +  {zgW-  zbB)  -  Ztp(W  -  B) 

where 

B  =  buoyancy,  Ibf 
W  =  weight,  in  air,  Ibf 
u  =  component  of  velocity  along  body  x  axis 
T'z  =  towcable  tension  along  body  x  and  z  axis 
systems 

Q  =  axial  force  coefficient  at  a  =  0 
Cjq  -  normal  force  coefficient  at  a  =  0 
=  pitching  moment  coefficient  at  a  =  0 
Qva  =  slope  of  the  normal  force  coefficient  at  a  0 
Cfna  =  slope  of  the  pitching  moment  coefficient  at 
a  =  0 

^  z  coordinates  of  center  of  buoyancy 
(C.B.),  ft 

Xg,  Zg=  X  and  z  coordinates  of  center  of  gravity  (C.G.),  ft 
xpp,  Ztp-  X  and  z  coordinates  of  the  towpoint  attachment,  ft 
a  =  angle  of  attack,  deg 
atrim  =  trim  angle  of  attack,  deg 

p  =  density  of  seawater,  1.9905  slug/ft^ 

A  =  reference  area  (86.59  ft^) 

/=  reference  length  (10.5  ft) 

4.0  Results  and  Discussion 

For  the  initial  VS  AERO  calculations,  the  TULP  ge¬ 
ometry  was  modeled  with  a  bluff  base.  The  launch  tube 
fairing  consisted  of  a  symmetric  airfoil.  The  approach  flow 
velocity  was  5  knots  (8.44  ft/s)  and  the  Reynolds  number 
(Rg )  was  7  million  based  on  a  reference  length  of  10.5  ft 
(launch  tube  diameter).  The  flow  angularity  was  0°.  Fig¬ 
ure  3  shows  an  isometric  rear  view  of  the  bluff  base  con¬ 
figuration.  Imposed  on  the  configuration  surface  are  the 
surface  pressure  coefficient  distribution  (Cp),  the  surface 
streamlines,  and  VSAERO  geometry  panels.  The  dark 
shaded  regions  represent  high-  and  low-pressure  regions. 
The  wake  shedding  off  the  bluff  base  is  represented  by 
the  panels  shown  to  the  left  of  and  aft  of  the  body.  The 
surface  streamlines  are  shown  as  streaks  along  the  sur¬ 
face  of  the  body.  The  surface  streamlines  terminate  at  loca¬ 
tions  where  VSAERO  determined  the  boundary-layer 
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shape  factor  to  be  of  a  value  that  indicates  flow  separa¬ 
tion.  Therefore,  the  regions  where  the  surface  streamlines 
have  terminated  approximately  coincide  with  regions  of 
flow  separation.  The  initial  bluff  base  calculations  revealed 
flow  separation  along  the  launch  tube  fairing  trailing  edge, 
fairwater  deck  trailing  edge,  and  sharp  leading  edge  of 
the  bow.  The  calculations  also  revealed  a  large  unsteady 
and  turbulent  wake  shedding  off  the  bluff  base. 

The  TULP  geometry  was  modified  to  alleviate  some 
of  the  flow  problems.  The  stern  of  the  main  deck  was 
lengthened  by  6  ft,  the  overall  structure  was  streamlined 
(ogive  shape)  over  the  aft  22  ft  of  the  main  deck,  and  the 
bow  was  rounded  along  the  leading  edge.  Additional 
streamlining  of  the  launch  tube  was  not  performed  at  this 
time.  The  modified  TULP  geometry  and  VSAERO  pan¬ 
eling  are  shown  in  Figure  4.  VSAERO  calculations  were 
then  performed  with  the  modified  geometry.  The  modi¬ 
fied  TULP  configuration  consisted  of  approximately  3500 
body  panels  and  770  wake  panels.  The  solutions  were 
obtained  by  using  a  14-block  formulation  for  the  Gauss- 
Seidel  iterative  solver.  The  solutions  were  coupled  with 
an  integral  boundary-layer  method  to  obtain  163  surface 
streamlines  and  boundary-layer  calculations.  Each  com¬ 
plete  VSAERO  calculation  included  6  wake  iterations  or 
relaxations  and  20  inviscid/viscous  iterations.  The  CPU 
time  for  each  solution  ranged  between  510  and  525  s. 

A  matrix  of  flow  conditions  for  a  steady  tow  speed 
of  5  knots  for  the  modified  TULP  was  defined  and  the 
computations  were  performed.  The  matrix  is  for  angles 
of  attack  (a)  of  0.0°,  ±1 .0°,  ±2.5°,  and  ±5.0°  each  at  angles 
of  sideslip  O)  of  0.0°,  ±1.0°,  ±2.5°,  and  ±5.0°.  Table  2 
lists  the  axial  (C^),  normal  (Q),  and  side  {Cy)  force  coef¬ 
ficients  as  well  as  the  pitching  moment  (C^),  rolling  mo¬ 
ment  (Q),  and  yawing  moment  (C„)  coefficients.  The  non¬ 
zero  sideslip  angle  computations  are  not  presented  in  this 
paper. 

Figure  5  shows  results  from  VSAERO  at  an  angle  of 
attack  of  +5°  for  the  modified  TULP  configuration.  The 
figure  shows  the  surface  streamlines  and  wake.  Compari¬ 
son  of  the  wake  of  the  modified  TULP  configuration  (Fig¬ 
ure  5)  and  the  bluff  body  TULP  configuration  (Figure  3) 

Table  2.  Hydrodynamic  force  and  moment  coefficients. 


shows  that  the  wake  off  the  streamlined  TULP  is  not  as 
unsteady  and  turbulent  as  the  wake  off  the  bluff  base 
TULP.  The  streamlining  of  the  stem  did  not  appear  to 
eliminate  flow  separation  along  the  trailing  edge.  How¬ 
ever,  additional  streamlining  may  help  reduce  or  elimi¬ 
nate  this  region  of  flow  separation.  Flow  separation  along 
the  bow  leading  edge  appeared  to  be  eliminated.  Addi¬ 
tional  streamlining  is  recommended  to  optimize  the  TULP 
configuration  for  reducing  flow  separation,  reducing  drag, 
and  further  reducing  the  turbulent  wake.  This  additional 
design  work  was  beyond  the  scope  of  this  study. 

Figures  6  through  9  depict  the  evolution  of  the  hy¬ 
drodynamic  forces  (C^,  C//,  Cy)  and  pitching  moment  (C^) 
as  a  function  of  the  number  of  VSAERO  iterations  for  the 
zero  angle  of  attack  and  zero  angle  of  sideslip  case.  In 
each  figure,  the  first  seven  data  points  are  the  inviscid 
calculations  that  include  one  rigid  wake  and  six  wake  re¬ 
laxations.  The  next  19  data  points  are  the  remaining  in¬ 
viscid/viscous  calculations.  These  remaining  iterations  are 
where  VSAERO  couples  the  integral  boundary-layer 
method  with  the  inviscid  calculations.  Oscillations  in  the 
coefficients  for  the  viscous  calculations  are  due  to  the 
application  of  the  transpiration  term  in  the  solution.  The 
transpiration  term  is  applied  in  the  boundary  condition  to 
account  for  the  boundary-layer  displacement  thickness. 
The  amount  of  transpiration  is  adjusted  on  each  iteration 
and  causes  small  differences  in  the  pressure  distribution 
from  iteration  to  iteration.  These  differences  in  the  pres¬ 
sure  distribution  show  up  in  the  integrated  hydrodynamic 
forces  and  pitching  moment.  These  oscillations  are  purely 
numerical  in  nature  and  are  not  an  expected  physical  phe¬ 
nomenon.  In  all  cases,  the  final  coefficient  value  is  taken 
at  the  last  inviscid/viscous  iteration. 

Figure  6  shows  the  axial  force  coefficient  (Q)  as  a 
function  of  iteration  number.  The  viscous  iterations  show 
a  large  increase  in  Q  due  to  addition  of  the  skin  friction 
and  transpiration  effects.  The  normal  force  coefficient  (Q) 
(see  Figure  7)  and  the  side  force  coefficient  (Cy)  (see 
Figure  8)  have  a  near-constant  inviscid  value.  The  vis¬ 
cous  iterations  show  a  small  decrease  in  the  magnitude  of 
the  normal  force  coefficient  and  a  small  increase  in  the 


a  (deg) 

Ca 

Cn 

c„ 

Ck 

Cl 

C„ 

-5.0 

0.3136 

-1.2040 

-7.0830 

0.0692 

0.2568 

-0.3968 

-2.5 

0.4105 

0.0515 

-4.3570 

0.0376 

0.2457 

-0.1623 

-1.0 

0.4586 

0.6288 

-2.3330 

0.0298 

0.3268 

-0.0383 

0.0 

0.4768 

1.0170 

-1.0060 

0.0430 

0.3430 

0.0467 

1.0 

0.4966 

1.5180 

0.1640 

0.0320 

0.2757 

0.1386 

2.5 

0.4897 

2.0020 

2.3640 

0.0197 

0.3833 

0.2741 

5.0 

0.4817 

2.9050 

5.8400 

0.0088 

02m 

0.5270 

Length  =  66  ft;  launch  tube  diameter 
C5a,=  10.5  ft;  5„f=  86.59  = 

=  10.5  ft;  velocity  =  5  knots  = 
10.5  ft;  a  =  ±5.0°;  p  =  0°. 

8.44  ft/s;  Reynolds  number  (R^ )  = 

7  million; 
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magnitude  of  the  side  force  coefficient.  The  small  increase 
in  side  force  may  be  due  to  an  effective  increase  in  asym¬ 
metry  due  to  the  displacement  thickness.  The  pitching 
moment  coefficient  (C^)  variation  with  iteration  number 
is  illustrated  in  Figure  9.  The  moment  reference  center  for 
these  calculations  was  chosen  to  be  along  the  lower  sur¬ 
face  of  the  main  deck,  near  the  juncture  point  of  the 
centerline  of  the  launch  tube  and  the  lower  surface  of  the 
main  deck.  The  inviscid  iterations  show  a  nearly  constant 
negative  value  of  the  pitching  moment  coefficient.  The  vis¬ 
cous  iterations  show  a  less  negative  (decrease  in  magni¬ 
tude)  pitching  moment  coefficient  due  to  a  smaller  viscous 
normal  force  coefficient. 

Figures  10  through  12  show  the  variation  of  the  hy¬ 
drodynamic  forces  as  a  function  of  angle  of  attack.  Each 
figure  shows  the  VS  AERO  results  at  zero  angle  of  side¬ 
slip  and  a  Reynolds  number  of  7  million.  Q  increases 
with  angle  of  attack  up  to  1°  and  then  decreases  with  a 
further  increase  in  angle  of  attack  (Figure  10).  C;v  increases 
nearly  linearly  between  ±5°  angle  of  attack  (Figure  11). 
There  appear  to  be  very  small  nonlinear  effects  in  the 
normal  force  coefficient.  The  asymmetric  shape,  which 
is  due  to  the  combination  of  the  fairwater  and  main  decks, 
is  most  likely  the  cause  for  the  small  nonlinear  effects  as 
a  function  of  angle  of  attack.  The  magnitude  of  the  side 
force  coefficient  (Cy)  is  non-zero  because  of  the  asym¬ 
metric  shape  of  the  TULP  configuration  (Figure  12).  Cy 
decreases  in  magnitude  with  increasing  angle  of  attack 
because  the  geometry  is  asymmetric  about  the  x-z  plane. 
The  asymmetry  is  due  to  the  lateral  offset  of  the  fairwater 
deck.  At  negative  angles  of  attack,  the  asymmetry  is  fully 
exposed  to  the  onset  flow,  resulting  in  a  larger  asymmet¬ 
ric  side  force.  However,  at  positive  angles  of  attack,  the 
asymmetry  is  partially  shielded  from  the  onset  flow,  re¬ 
sulting  in  a  smaller  asymmetric  side  force. 

The  variation  of  the  hydrodynamic  moments  with 
angle  of  attack  is  portrayed  in  Figures  13-15.  Once  again, 
the  VS  AERO  results  are  for  zero  angle  of  sideslip  and  a 
Reynolds  number  of  7  million.  Figure  13  shows  that 
increases  linearly  between  ±5°  angle  of  attack.  There  are 
small  nonlinear  effects  near  ±5°  angle  of  attack.  The  non¬ 
linear  effects  arise  from  the  nonlinear  effects  of  normal 
force,  which  is  likely  due  to  the  asymmetric  shape  of  the 
modified  TULP  configuration.  The  positive  slope  of  the 
pitching  moment  suggests  an  unstable  configuration. 
However,  when  buoyancy,  weight,  and  tow  forces  are  con¬ 
sidered,  the  TULP  is  stable,  as  shown  in  a  later  section  of 
this  paper.  The  rolling  moment  coefficient  (Q)  is  non¬ 
zero  at  all  angles  of  attack  and  varies  between  0.25  and 
0.39  (Figure  14).  Q  is  non-zero  because  of  the  geometric 
asymmetry  about  the  x-z  plane.  The  yawing  moment  co¬ 
efficient  (C„)  increases  linearly  between  ±5°  angle  of  at¬ 
tack  (Figure  15). 

The  VS  AERO  computed  coefficients  were  used  in 


conjunction  with  mass,  buoyancy,  and  towpoint  estimates 
to  compute  the  trim  angle  of  attack  (atrim)  for  the  refined 
TULP  configuration,  under  steady-state  tow.  The  trim 
angle  of  attack  was  computed  for  three  different  states 
during  the  underwater  launch  sequence  of  a  Trident-size 
Launch  Test  Vehicle  (LTV).  The  three  states  are  LTV  in 
launch  tube  (prior  to  launch),  launch  tube  empty  (just  af¬ 
ter  launch),  and  the  launch  tube  flooded  with  seawater 
(after  launch).  The  empty  launch  tube  represents  the 
“worst  case”  in  terms  of  platform  stability  due  to  the  re¬ 
duced  submerged  metacentric  height  (i.e.,  the  distance 
between  the  center  of  gravity  and  the  center  of  buoyancy). 
Under  all  three  conditions,  the  calculations  show  that  the 
modified  TULP  configuration  under  steady  tow  will  trim 
at  an  angle  of  attack  less  than  or  about  1.0°  (Table  3). 


Table  3.  Platform  stability  under  steady  tow. 


Submerged 

Trim  Angle 

Metacentric 

of  Attack 

Launch  State 

Height  (ft) 

(deg) 

LTV  in  launch  tube 

4.6 

1.03 

Launch  tube  empty 

1.3 

0.66 

Launch  tube  flooded 

4.3 

0.97 

The  stability  of  the  modified  TULP  is  due  to  the  place¬ 
ment  of  the  buoyancy  pontoons  along  the  deck  and  the 
lead  ballast  at  the  bottom  of  the  launch  tube.  While  this 
study  shows  stable  steady-state  tow  performance  with 
small  angles  of  attack,  a  more  comprehensive  dynamic 
model  would  need  to  be  developed  to  explore  the  full  range 
of  platform  tow  motions. 

Another  important  aspect  of  the  TULP  concept  is  its 
time  history  response  to  an  LTV  test  launch.  Launching 
an  LTV  or  missile  involves  a  gas  or  steam  generator  with 
an  impulse  load  transmitted  to  the  launching  platform. 
The  platform  response  must  be  sufficiently  small  to  en¬ 
sure  that  the  platform  recoil  does  not  alter  the  trajectory 
of  the  LTV  or  corrupt  the  collection  of  vital  engineering 
data.  The  recoil  of  the  platform  is  a  function  of  the  launch 
impulse  as  well  as  the  mass  and  added  mass  of  the  plat¬ 
form.  The  TULP  has  a  large  hydrodynamic  added  mass 
in  heave  due  to  the  shape  and  size  of  the  main  deck. 

In  the  concept  study  of  Ref.  1 ,  a  transfer  function 
approach  was  used  to  estimate  the  platform  recoil  during 
a  launch  of  a  Trident-size  LTV.  The  analysis  showed  that 
the  platform  would  initially  recoil  downward  approxi¬ 
mately  0.65  ft  with  a  vertical  velocity  of  about  1.4  ft/s.  In 
comparison,  the  LTV  will  have  traveled  an  in-tube  dis¬ 
tance  of  40  ft  and  will  have  attained  a  velocity  of  nearly 
60  ft/s.  The  response  of  a  Trident  submarine  to  a  launch 
will  be  a  small  fraction  of  that  of  the  test  platform. 
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However,  the  difference  in  their  respective  response  char¬ 
acteristics  should  not  be  significant  in  terms  of  assessing 
the  in-tube  and  underwater  flight  performance  of  the  LTV. 

Subsequent  to  the  initial  launch  impulse,  the  plat¬ 
form  will  experience  a  light  condition  (launch  tube  empty) 
followed  by  a  gradual  return  to  a  constant  weight  trim  as 
seawater  floods  the  launch  tube.  A  detailed  transient  analy¬ 
sis  with  a  6  degree-of-freedom  simulation  would  be  nec¬ 
essary  to  fully  understand  all  of  the  response  characteris¬ 
tics  to  a  launch  impulse.  The  USAERO  code,®  a 
time-dependent  panel  method  code,  can  be  used  to  evalu¬ 
ate  the  TULP  response  to  a  launch  impulse.  USAERO 
solves  for  the  inviscid  potential  flow  for  single  or  mul¬ 
tiple  bodies  in  general  motion.  The  USAERO  code  is  a 
sister  code  to  the  VSAERO  code  and  was  initially  devel¬ 
oped  from  an  early  version  of  VSAERO.  The  VSAERO 
geometric  inputs  are  the  same  as  the  geometric  inputs  re¬ 
quired  for  USAERO. 

5.0  Conclusions 

A  point  design  for  a  TULP  is  reported  in  this  paper. 
The  TULP  concept  has  been  configured  to  provide  a  test 
bed  that  replicates  the  hydrostatic  pressure  (depth)  and 
the  submarine  deck  flowfield  (crossflow)  at  full  scale  to 
capture  the  complex  phenomenology  of  underwater 
launch.  The  configuration  provides  a  large  submerged 
metacentric  height  for  pitch  and  roll  stability,  is  stream¬ 
lined  in  the  fore-aft  direction  for  low  forward-way  drag, 
and  has  a  high  added  mass  in  heave  to  control  launch  re¬ 
coil  motions. 

VSAERO  panel  code  calculations  were  performed 
on  the  initial  TULP  concept  geometry  to  assess  flow  sta¬ 
bility  and  hydrodynamic  performance.  These  calculations 
revealed  regions  of  flow  instability  and  separation.  The 
TULP  geometry  was  revised  to  correct  the  flow  problems, 
and  then  the  complete  hydrodynamic  forces  and  moment 
coefficients  were  computed.  The  hydrodynamic  coeffi¬ 
cients  were  combined  with  information  about  mass,  buoy¬ 
ancy,  and  towpoint  location  to  show  a  small  trim  angle  of 
attack  under  a  variety  of  launch  tube  conditions.  These 
results  are  encouraging  indicators  that  the  TULP  concept 
could  provide  a  viable  test  platform. 

6.0  Recommendations  for  Future  Work 

The  recoil  of  the  platform  to  a  launch  impulse  can  be 
addressed  with  the  USAERO  code,  an  unsteady  panel 
method  code.  USAERO  uses  the  same  geometry  inputs 
as  the  VSAERO  code  and  is  a  natural  follow-on  to  the 
VSAERO  results  reported  here.  This  paper  only  addressed 
the  hydrodynamic  stability  of  the  platform  in  the  pitch 


plane.  It  is  recommended  that  additional  work  be  per¬ 
formed  to  address  the  tow  stability  with  non-zero  angles 
of  sideslip.  Also,  a  6  degree-of-freedom  simulation  that 
includes  the  tow-cable  dynamics  should  be  performed. 
Model-scale  tests  of  the  TULP  in  a  tow  tank  should  be 
carried  out  to  measure  the  forces  and  moments  acting  on 
the  platform. 
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Axial  Force  Coefficient 


Iteration 


Figure  6.  Axial  force  coefficient  as  a  function  of  iteration 
number. 


Figure  3.  Launch  platform  (before  geometry  modifica¬ 
tions). 


Figure  4.  Modified  TULP. 


Figure  5.  Boundary  layer  and  wake  flow  for  the  modified 
TULP, 


Norma!  Force  Coefficient 
«  =  0.0  deg..  Re  =  7  Million  CN| 


1.50  . 

1.25  - 

1  * 

i  ^ 

0.75  - 

Inviscid  - 

^  1  ►  ' 

Viscous  _ 

0.50  - 

1  1 

0.25  - 

0.00  - 

0  S  10  15  20  25  30 


Iteration 

Figure  7.  Normal  force  coefficient  as  a  function  of  itera¬ 
tion  number. 
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Figure  8.  Side  force  coefficient  as  a  function  of  iteration 
number. 
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Figure  9.  Pitching  moment  coefficient  as  a  function  of 
iteration  number. 
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^  ^  ^  .  o  Figure  13.  Pitching  moment  coefficient  as  a  function  of  a. 

Figure  10.  Axial  force  coefficient  as  a  function  of  a. 
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Figure  1 1 .  Normal  force  coefficient  as  a  function  of  a. 
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Figure  14.  Rolling  moment  coefficient  as  a  function  of  a. 
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Figure  12.  Side  force  coefficient  as  a  function  of  a.  Figure  15.  Yawing  moment  coefficient  as  a  function  of  a. 


9 

American  Institute  of  Aeronautics  and  Astronautics 


DEPARTMENT  OF  THE  NAVY 

STRATEGJC  SYSTEMS  PROGRAMS 
1931  JEFFERSON  DAVIS  HWY 
ARLINGTON.  VA  22241*5362 


9199 


5531 

SP1622 

Ser  U091896027 


1.0 


8rP 


From:  Director,  Strategic  Systems  Programs 
To:  Director,  The  Johns  Hopkins  University 

Applied  Physics  Laboratory 
(Attn:  L.  Croston) 

Johns  Hopkins  Road 
Laurel,  MD  20723-6099 

Sub j :  CLEARANCE  OF  MATERIAL  FOR  PUBLIC  RELEASE;  SSP  SECURITY 
REVIEW  CASE  96-55 


Ref:  (a)  JHU/APL  Itr  96-173  of  29  Aug  96 

End:  (1)  Paper,  "Development  of  a  Towed  Underwater  Launch 

Platform" 


1.  Enclosure  (1),  forwarded  by  reference  (a),  has  been  reviewed 
and  is  approved  for  unlimited  distribution. 


P.  M.  WALLACE 
By  direction 


Copy  to:  (w/o  end) 

LMMS,  Sunnyvale,  (T.  Tucci) 


xc:  Parthasarathy 

L •  Levy 
J.  Gibson 


Im 


091896027/000 


strategic  systems  programs 

WASHINGTON.  DC  20376-5002 


DUG  0UAI.<TTY  IFc[^EC‘ClD  I 


CJD 


